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Abstract 

We construct generalizations of the Calogero-Sutherland-Moser system by 
appropriately reducing a model involving many unitary matrices. The result- 
ing systems consist of particles on the circle with internal degrees of freedom, 
coupled through modifications of the inverse-square potential. The coupling 
involves SU{M) non-invariant (anti)ferromagnetic interactions of the internal 
degrees of freedom. The systems are shown to be integrable and the spectrum 
and wavefunctions of the quantum version are derived. 
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The inverse- square interacting particle system []I|, |^ and its spin generaliza- 
tions 0, |, H, 13, H, H are important models of many-body systems, due to their 
exact solvability and intimate connection to spin chain systems |TT], |r^, |TB 



2-dimensional Yang-Mills theories |14, 16 1 etc. 



A useful route for studying these systems is through reductions of appropriate 
matrix models fl^, |18|]. Although this does not give rise to all known models in 



the quantum case, it is still a very direct way to show integrability and solvabil- 
ity. In this paper, we will explore this approach further and show how, starting 
with many matrices effectively coupled through constraints, we can derive further 
generalizations of these systems with internal degrees of freedom. 
The starting point will be the many-matrix lagrangian 

M 1 

The Ui are N x N unitary matrices depending on time t and overdot stands for time 
derivative. The at are real positive parameters. For convenience, we shall normalize 
the Qi such that 

M 

E«^ = l (2) 

The above system is clearly integrable and solvable, being nothing more than 
a collection of independent unitary matrix models. In the absence of further con- 
straints it would give rise to a collection of independent spin-generalizations of the 
Sutherland model. It is, however, possible to reduce it in a different way: the above 
model is invariant under (independent) left- and right-multiplications of the Ui by 
constant matrices. This manifests in the existence of conserved quantities, namely 

Ri = --Ur'ili , L, = -U.Ur^ (3) 
tti ai 

The above conserved hermitian matrices Poisson-commute with themselves to the 
U (N) algebra and mutually to zero. We will choose to fix the values of the following 
time-independent quantities: 

- -U-'if, + -^IJ,^,Ur^\ = P, (4) 

where we have adopted periodic conditions in i, that is, M + 1 = 1, = M. Clearly 
the chosen expressions for the Pi, being the sum of independent U{N) generators, 
also constitute mutually commuting Poisson-?7(A^) matrices. 

The key to a successful reduction of this system to particles lies in an appropriate 
parametrization of the Ui. We will choose 

Ui = Wi'_\A^^W, (5) 
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where Wi are unitary matrices and A = diag{e'^^") is a diagonal unitary matrix. The 
reason for the choice of exponents in (Q) will be apparent in the sequel. From the 
relation 

U = Ui---Un = Wj^^AWn (6) 

we see that Wn and e*^" are the diagonalizer and the eigenvalues of the matrix U, and 
the remaining Wi are determined recursively from (^). The above parametrization 
has the redundancy generated by left-multiplication of all Wi by the same diagonal 
matrix. This will lead to a 'gauge constraint' later on. 

With the above parametrization the expressions for the lagrangian and the Pi 
become, after some algebra, 

^ = t^E ^ |-(A-^A)^ +(- + - -A_iA'^«AA-'^4 (7) 

WiPiWr^ =Ki=[ — ^ d - -A-"»£,_iA"> - ^A'^'+i£i+iA-'^>+^ (8) 

where we defined 

WiWr' = iL, (9) 

From the structure of the above lagrangian, we see that the generator of left-rotations 
for Wi (found by differentiating L with respect to iWiWj~^) is exactly Ki as expressed 
in (P). Therefore, the Ki are also mutually commuting Poisson-f/(A^) generators. 
This justifies the choice of exponents in (|^), since any other choice would spoil the 
decoupling and U{N) nature of the Ki. We stress, however, that the Ki are no more 
conserved and are, thus, dynamical quantities. Further, from (|D we see that the 
diagonal elements of Ki must satisfy 

M 

Kinn = for each n (10) 

1=1 

which is the 'Gauss law' originating from the 'gauge invariance' of the parametriza- 
tion (^ as stated previously. The off-diagonal elements of satisfy 

M 

^ ^ -M^tj (^^mn)-^7'mn -^imn (-^-^) 

i=i 

where 6mn = dn and the matrix M(x) is defined as 

M,,(x) = (- + 6i, - -e-'^»^'5,,,+i - ^e'^'+i"5,+i,,- (12) 

\ai aj_|_i J tti aj_|_i 

(again, taking 5j,Af+i = In terms of the matrices Ki and the angles 6'„ the 

hamiltonian (which is identical to the lagrangian since it only contains kinetic terms) 
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becomes 



n=l m^n=l « j'=l n=l j j'=l 

(13) 

where p„ = 6'^ is the canonical momentum of 6n and M(0) is the matrix M(0) 
projected to the subspace orthogonal to its zero-eigenvalue eigenvector. 

In order to invert the matrix M, we note that it is essentially a discrete second 
derivative on a set of points at distances a^, conjugated with an exponential factor. 
Generalizing from the case of equally-spaced points (that is, the Cartan matrix of 
SU{M)), where the inverse is proportional to the distance |i — j|, we define 



and try a form 



= E ' ^ij = h- bj (14) 



Mr^ix) = e-''^" {A{x) + B{x)h, + C(x)|%|} e''^" (15) 



The elements with z, j either 1 or M need special attention due to the 'wrap around' 
properties of M. (Note that hi is not periodic in its index, but rather 6i+Af = h + l.) 
This fixes the coefficients A,B,C and the inverse of M becomes 



The matrix M(0) has a zero eigenvalue corresponding to the eigenvector (1, . . . 1). 



In order to compute the inverse of M it suffices to take the limit a; ^ in (|16D 
and discard pieces that are constant in i,j since they clearly project on the zero- 
eigenvalue eigenspace. We find 



mxr% = -^2 + ^ + + (17) 



Discarding the constant piece and taking x ^ we obtain M ^ . Putting everything 
together, we finally obtain 

with the potentials Vij and Vij defined as 

- e ]^^^.^, I + ' 2 '° 2 2 J ' 2 ^^^^ 

This is a generalization of the Sutherland model of particles on the circle with spin 
degrees of freedom, encoded by the Ki, where the different Ki couple through a 
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parametric modification of tlie inverse-sine-squared potential. Interestingly, this 
is the same classical model as the one obtained by Blom and Langmann starting 
from the 2-dimensional Yang-Mills point of view JTUj. The standard spin-Sutherland 
model is recovered either as the limit where all Ki but one are zero, or as the limit 
where all a, but one are zero, in which case only the sum of all Ki appears in the 
hamiltonian. 

The quantum version of this model proceeds along similar lines. The constant 
matrices Pi become, now, U{N) generators transforming under some fixed repre- 
sentations of SU{N) and carrying some U{1) charge (since they are not necessarily 
traceless). We shall choose them to be irreducible, the reducible case being simply 
the direct sum of models with one irreducible component for each Pi. Therefore, the 
model is labeled by a set of irreps and charges qi carried by the Pi. The matrices 
Ki become time- varying SU (N) generators in the same irrep and with the same 
charge qi as Pi. In a more standard notation, putting 

K,= Y: T'^Kt + f (20) 

a=l 

with T"" the fundamental SU{N) generators, the Ki obey the SU{N) algebra for 



each i and commute for different i, while the qi are central. The hamiltonian (|T 
remains valid for the quantum operators as well; the only term requiring ordering, 
namely VnKimnKinmi is automatically symmetrized by the summation over m, n. 
To turn the Ki into internal degrees of freedom we follow the standard construc- 



tion [|I^ (see also ^^). The Ki can be realized a la Jordan- Wigner in terms of 
bosonic oscillators 

1 qi 

Kimn = {Mkm.^ikn — -TJ^mn ^Ikl^ikl) + T7^mn (21) 
k=l 1=1 

where di is the number of rows in the Young tableau of r^. The Aikm are a collection 
of commuting bosonic ladder operators 

[-Aikmy ^lln] ~ ^ij^kl^mn (22) 

The above imbeds r^, as well as all other irreps with up to di rows, in the Fock 
space of Aikm- To simplify the notation and interpretation of the model we will 
choose di = 1 for all K^. We can always achieve the case di > 1 as a model with 
M' = Y^iti di matrices with all di = 1. Choosing 

i—l i 

a,=0 for Y.dk<J<Y.dk (23) 

k=l k=l 

makes all the Kj in this range, overall di in number, appear only through their 
sum in the hamiltonian, reproducing the original di > 1 matrix. We thus drop the 
summation and index k in (^) and (^2]) . 
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The gauge constraint (|iy) in terms of (^Tj) implies 

M M 

A\^Ain = constant = m for all n, X! ~ ^ (^^) 

1=1 i=i 

From the above we see that the Fock states generated by the M oscillators Ai^ for 
each n transform in the m-fold symmetric irrep of SU{M). We therefore define the 
N mutually commuting SU{M) generators 

Snij = Aj^Ajn — J^^ij (25) 



carrying the m-fold symmetric irrep of SU{M). In terms of these, substituting ( pi]) 
in ([T^), we obtain after some algebra 

H = - ^pI + -Yl iYVij{d„,n)SmijSnji + Vo{6mn) ^ ^ 

" m^n \ ij 

+ 2 ^ ^ ^jj'S'njjS'njj + ^ Vij {qiQj — SiiSjj) (26) 

n ij ij 

where we defined the total spin S = Y.n=i and Vo{x) = Vii{x) = 1/(2 sin |)^. The 
above hamiltonian refers to a system of particles with m-fold symmetric SU{M) 
spins 5"^. Due, however, to the existence of the matrix V{x) in the coupling of spins, 
the interactions above are not S'f/(M)-invariant (the total spin S is not conserved). 
Only the diagonal elements (that is, the Cartan part) of the total spin are conserved. 
The appearance of the total spin in (^) can be eliminated by choosing the charges 



mA^ , , 

Qi = Sii = ni- -j^ (27) 

where n, is the number of boxes in r^. (In this construction the rj are symmetric, 
so Uj are their lengths.) This is a natural choice, identifying the U{1) charge to the 

charge of rj and subtracting the average Zj^ charge of all irreps. 

In terms of the original matrix problem, the hamiltonian is the sum of indepen- 
dent laplacians over the space U{N). The U{1) part will lead to charges Qi (the 
momentum of the 'center-of-mass' of each Ui) which, by (|^) and (^), have to satisfy 

Qi - Qi+i =qi , J2 ^iQi = J2Pn = P (28) 

i n 

This fixes the Qi in terms of the qi and the total momentum P. Factoring out 
the U{1) part for each Ui, we are left with laplacians on the SU{N) manifold. It 
is known that the matrix elements of each irrep R of SU{N) form a degenerate 
eigenstate multiplet of the laplacian with eigenvalue given by the quadratic Casimir 
C2{R) of the irrep. Therefore, the eigenstates of the matrix hamiltonian are 

Q 

-^({Ui}) = n(detf/^)^i?^(f/^)a,ft (29) 
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where Ri are irreps of SU{N) {ai, f3i label their matrix elements). The energy 
eigenvalue E and degeneracy D corresponding to this state are 



E = J:^{C,{R,) + $\ , D = Y[{d^mR,r (30) 




N 

J l—L 

Each Ri{Ui)aii3i transforms under left-rotations of Ui in the Ri irrep (acting on the 
left index) and under right-rotations of Ui in the conjugate irrep Ri (acting on the 
right index). The condition (^), however, tells us that the sum of the right-generator 
for Ui and the left-generator for [/j+i are constrained to be in the irrep of SU{N) 
carried by Pj. Therefore, we must project the states Pi and Oj+i, transforming in 
the Ri and Ri+i respectively, to the subspace of states transforming in the rj. Call 
G(Pj, aj+i|rj, 7j) the Clebsch-Gordan coefficient that projects these states 
to the 7i state of r^. Then the final energy eigenstate wavefunction for this model 
becomes 

^ Q- 

^^i{U^})= E n(det?7,)^P.(f/.)a,ftG'(i?„A;i?^+i,«mki,7^) (31) 

The indices 7, do not imply a degeneracy of states corresponding to the dimension- 
ality of Ti, since these are states fixed by the constraints (^) and are not summed 
over. The degeneracy of each state labeled by the Ri is given by the number of 
times that the irreps are contained in the direct products Ri x Ri+i or, equiva- 
lently, the number of times that Ri+i is contained in Ri x r^. Calling this integer 
D{Ri, Vi] Ri+i), we have for the total degeneracy 

M 

D = Y[D{R,,r,;Ri+,) (32) 

i=l 

Clearly, an irrep Ri can appear in \E' only if it is contained in the product Ri x 
ri ■ ■ ■ X r^. This requires that the total Z^v charge of the r's be zero, that is, 
the total number of boxes in the Young tableaux of rj should be a multiple of N. 
This is indeed the case, the total number being n = mA^. The above, upon proper 
reinterpretation, gives the spectrum and wavefunctions in each corresponding sector 
of the particle-spin model as classified by the Cartan elements of the total spin. 

We conclude by mentioning that a construction similar to (^TJ) in terms of 
fermionic oscillators would give rise to a model as in ( p6D but with negative sign 
(ferromagnetic) spin interactions. 

I would like to thank E. Langmann for comments and for communicating his 
results prior to publication. 
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